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Abstract 

Suppose we have an n-qubit system, and we are given a collection of local density 
fT^ ' matrices pi, . . . ,pm, where each pi describes a subset d of the qubits. We say that 

^ I the Pi are "consistent" if there exists some global state a (on all n qubits) that 

J^ ' matches each of the pi on the subsets Ci. This generalizes the classical notion of 

the consistency of marginal probability distributions. 

We show that deciding the consistency of local density matrices is QMA-complete 
(^ [ (where QMA is the quantum analogue of NP). This gives an interesting example of a 

^^ • hard problem in QMA. Our proof is somewhat unusual: we give a Turing reduction 

^P I from Local Hamiltonian, using a convex optimization algorithm by Bertsimas and 

Vempala, which is based on random sampling. Unlike in the classical case, simple 
^1^1 mapping reductions do not seem to work here. 

-(— > 

c, 

^ , 1 Introduction 

a'' 

Quantum mechanical systems exhibit many unusual phenomena, such as coherent super- 
positions and nonlocal entanglement. It is interesting to compare this with the behavior 
p\ ' of classical probabilistic systems. In a classical system, such as a Markov chain or a 

j^ ■ graphical model, one may have correlations or dependencies among different parts of 

the system; in particular, local properties can affect the joint probability distribution of 
the entire system. Many quantum systems have a similar flavor, though their behavior 
is more complicated. In this paper, we investigate one problem of this kind, and its 
relationship to the complexity class QMA. 

First, consider a classical problem. Suppose we have random variables Xi, . . . ,Xn, 
with some unknown joint distribution D, and we are given marginal distributions Di, . . . , D„ 
where each Di describes a subset Ci of the variables. (We assume that the random vari- 
ables Xj take on values in some fixed finite set, and the subsets Ci have size at most 
some constant k.) Does there exist a joint distribution D that matches each of the 
marginals Di on the subsets Ci? If so, we say that the marginals Di are "consistent." 

Deciding the consistency of marginal distributions is NP-hard, by a straightforward 
reduction from 3-coloring. (We are given a graph G = {V,E). For each vertex v £ V, 
construct a random variable Xy which takes on values in {r, g, b}. For each edge {u, v) G 



E, specify that the marginal distribution of X^ and X^ must be uniform over the set 
{r,g, b}^ \ {rr,gg, bb}. These marginals are consistent iff G is 3-colorable.) 

Now consider the generalization of this problem to quantum states. (This problem 
was first suggested to me by Dorit Aharonov, in connection with the class QCMA [Ij.) 
Suppose we have an n-qubit system, and we are given local density matrices pi, . . . , p^, 
where each pi describes a subset Cj of the qubits. Does there exist a global state a on 
all n qubits that matches each of the local states pi on the subsets Cj? If so, we say that 
the local states pi are "consistent." 

We will show that this problem is QMA-complete, where QMA is the quantum 
analogue of NP. QMA is the class of languages that have poly-time quantum verifiers, 
where the witness is allowed to be a quantum state. QMA arises naturally in the study 
of quantum computation, and it also has a complete problem, Local Hamiltonian, which 
is a generalization of A:-S AT [21 [3] . 

Our result is interesting, because we only know of a few QMA-complete problems, 
and most of them look like universal models of quantum computation. For instance, 
the fact that Local Hamiltonian is QMA-complete [21 [U O [6] is closely related to the 
fact that adiabatic quantum computation is equivalent to the standard quantum circuit 
model |7]. Other QMA-complete problems such as Identity Check involve properties of 
quantum circuits [8]. The Consistency problem, however, does not seem to embody any 
particular model of quantum computation; this will become clearer when we present our 
reduction from Local Hamiltonian. 

Why are there so few QMA-complete problems, when there is such an astonishing 
variety of NP-complete problems? The reason seems to be that the techniques used to 
show NP-hardness, such as mapping reductions using combinatorial gadgets, break down 
when we apply them to a "quantum" problem like Local Hamiltonian. For instance, to 
reduce Local Hamiltonian to the Consistency problem, we would try to use local density 
matrices to "simulate" local Hamiltonians. But we run into problems due to the presence 
of non-commuting matrices. (In cases where quantum gadgets do work, such as [3 [6], 
they are much more subtle than classical gadgets.) 

Instead, our proof that the consistency problem is QMA-hard uses a randomized 
Turing reduction from Local Hamiltonian. The basic idea is that Local Hamiltonian can 
be expressed as a convex program in polynomially many variables, which can be solved 
using convex optimization algorithms, given an oracle for the Consistency problem. In 
particular, we use a class of convex optimization algorithms [9l ITOl [TT] which are based 
on random walks, and only require a membership oracle, rather than a separation oracle. 
We also use a nifty representation of local density matrices in terms of the expectation 
values of Pauli matrices. 

Note that the Consistency problem has a rather different structure from Local Hamil- 
tonian. For instance, a local density matrix contains complete information about the 
local state of the system, whereas in many cases a local Hamiltonian only constrains the 
local state of the system to lie within a certain subspace. 

Finally, we remark that our reduction from Local Hamiltonian to Consistency pre- 
serves the "neighborhood structure" of the problem, in that the local density matrices 
act on the same subsets of qubits as the local Hamiltonians. So, using the QMA-hardness 
results for 2-Local Hamiltonian ^ and Local Hamiltonian on a 2-D square lattice [6], 



we can immediately get QMA-hardness results for the corresponding special versions of 
the Consistency problem. 

We also mention some related work. In [13], one considers the Common Eigenspace 
Problem, verifying the consistency of a set of eigenvalue equations Hi\tp) = Xi\^p), where 
the operators H^ commute. We do something similar, translating each local density 
matrix into constraints on the expectation values of Pauli matrices, though in our case 
the Pauli matrices do not commute. Also, in [T^, one considers a quantum analogue of 
2-SAT, where we seek a state {ip) whose local density matrices have support on prescribed 
subspaces. However, this problem is more closely related to Local Hamiltonian than to 
Consistency, since the constraints can be written in the form HjI^) = where the Hi 
are local projectors. 

2 Preliminaries 

2.1 Density Matrices 

A quantum state of an n-qubit system is represented by a density matrix, which is a 
2" X 2" positive semidefinite matrix with trace 1. A classical joint probability distribution 
on n bits is a special case, where the density matrix is diagonal, and the diagonal entries 
are the probabilities of the 2" possible outcomes. A subset of qubits C is described by 
a reduced density matrix, which is obtained by taking the partial trace over the qubits 
not in C. This is analogous to a marginal distribution, which is obtained by summing 
over some of the variables. 

We measure the difference between two quantum states using the Li matrix norm, 
11/9 — (t||i = tr\p — a\. Note that this is also called the trace or statistical distance (when 
normalized by a factor of 1/2). 

Let X, Y and Z denote the Pauli matrices for a single qubit, and define V = 
{I,X,Y,Z}. We can construct n-qubit Pauli matrices by taking tensor products P = 
Pi (g) • • • (g) P„ E -p^n^ Any 2^-dimensional Hermitian matrix can be written as a real 
linear combination of n-qubit Pauli matrices. Furthermore, the n-qubit Pauli matrices 
are orthogonal with respect to the Hilbert-Schmidt inner product: tr{P^Q) = 2" if 
P = Q, and otherwise. So, if a is an n-qubit state, we can write it in the form 



l-n E «-^' 



2 

where the coefficients are uniquely determined by ap = tr (Pa); note that these are the 
expectation values of the Pauli matrices P. This application of the Pauli matrices is 
closely related to quantum state tomography. 

2.2 QMA and the Local Hamiltonian Problem 

The class QMA, or "Quantum Merlin- Arthur," is defined as follows [2^, ^: a language L 
is in QMA if there exists a poly-time quantum verifier V and a polynomial p such that 

• If X G L, then there exists a quantum state p on p(|x|) qubits such that V{x,p) 
accepts with probability > 2/3. 



• If X ^ L, then for all quantum states p on p(|x|) qubits, V{x,p) accepts with 
probability < 1/3. 

(Here, |x| denotes the length of the string x.) This is similar to the definition of NP, 
except that the witness is allowed to be a quantum state, and the verifier is a quantum 
circuit with bounded error probability. 

The Local Hamiltonian problem is defined as follows: 

Consider a system of n qubits. We are given a Hamiltonian H = Hi + • • • + 
Hm-, where each Hi acts on a subset of qubits Ci C {1,... ,n}. The Hi 
are Hermitian matrices, with eigenvalues in some fixed interval (for instance 
[0, 1]), and each matrix entry is specified with poly(n) bits of precision. Also, 
m < poly(n), and each subset Ci has size \Ci\ < k, for some constant k. 

In addition, we are given two real numbers a and b (specified with poly(n) 
bits of precision) such that b — a > l/poly(n). 

The problem is to distinguish between the following two cases: 

• li H has an eigenvalue that is < a, output "YES." 

• If all the eigenvalues of H are > 6, output "NO." 

Kitaev showed that Local Hamiltonian is in QMA, and the case of /c = 5 is QMA-hard 
[21 [3]. With greater effort, one can show that Local Hamiltonian with k = 2 is also 
QMA-hard [3l|5]. 

2.3 Convex Programming 

Consider the following version of convex programming: 

Let K C M" be a convex set, specified by a membership oracle Ox- 
Let f : i^ ^ M be a convex function, which is efficiently computable. 
Find some x £ K that minimizes f{x). 

Note that the membership oracle Ok is not as powerful as a separation oracle. We 
would like to solve this problem with precision e; that is, we want to find some x that 
lies within distance e of an optimal solution x*. 

We can solve this problem in time poly(n, log(l/e)), using an algorithm by Bertsimas 
and Vempala which is based on random sampling [Hill]. Actually, for our purposes we 
only need to solve the special case where / is a linear function; for this case, we can 
use a slightly faster simulated annealing algorithm [lOj, or an algorithm based on the 
shallow-cut Ellipsoid method [12j. But for simplicity we will stick with the Bertsimas 
and Vempala algorithm. 

Theorem 1 (Bertsimas and Vempala) Consider the convex program described above. 
Suppose K is contained in a ball of radius R centered at the origin. Also, suppose we 
are given a point y, such that the ball of radius r around y is contained in K . Then this 
problem can be solved in time poly(n,L), where L = \og{R/r). 



3 Consistency of Local Density Matrices 

We define the Consistency problem as follows [Ij: 

Consider a system of n qubits. We are given a collection of local density 
matrices pi, . . . , pm, where each pi acts on a subset of qubits Cj C {1, . . . , n}. 
Each matrix entry is specified with poly(n) bits of precision. Also, m < 
poly(n), and each subset Ci has size |Cj| < k, for some constant k. 

In addition, we are given a real number f3 (specified with poly(n) bits of 
precision) such that /3 > l/poly(n). 

The problem is to distinguish between the following two cases: 

• There exists an n-qubit state a such that, for all i, Ijtrn _„i_(^. (cr) — 
Pi 111 = 0. In this case, output "YES." 

• For all n-qubit states u, there exists some i such that \\ti\i,, ,,,n}-Cii'^) ~ 
Pi 111 ^ P- I^ this case, output "NO." 

Theorem 2 Consistency is in QMA. 

Proof: The basic idea is as follows. Given a witness state a, the verifier will pick a 
subset Ci, and perform measurements to compare a (on the subset d) to pi. There is 
a complication, however, because the verifier requires many independent copies of the 
witness a, and the prover might try to cheat using entanglement among the different 
copies. One can deal with this problem using a Markov argument. 

A concise way to present this result is to use QMA+, an alternate characterization 
of QMA [15j . We construct a QMA+ "super- verifier" for Consistency: 

Choose i € {1, . . . , m} at random. Choose a Pauli matrix Q € V^' ^' (acting 
on the subset Ci), at random. 

Output a quantum circuit M which does the following: measure the observ- 
able Q <Si I on the state a; if the measurement outcome is -|-1, accept; if it is 
— 1, reject. [So the circuit accepts with probability | + ^tr((Q (g) /)<t).] 

Set the target probability to be r = ^ + ^ tr(Qpj), and set the error threshold 
to be s = i(/3/4l^'l). [Note that s > |(/3/#) = l/poly(n).] 

To see why this works, write the difference between a and pi on subset C, as follows: 

tr{i,...,n}_c,(f^) -Pi = ^|c^ Yl ^Q*5, 

where 'jq = tr((Q 1)a) — iv{Qpi). On a "yes" instance, we always have 7q = 0, so the 
circuit M accepts with probability r. On a "no" instance, there must be some choice of 
Ci and Q such that |7q| > /?/4l'^*l; to see this, use the triangle inequality to write: 

/? < l|tr{i,...,n}-c,(cr) - PiWi < Yl I'^qI- 

In this case, the circuit M accepts with probability that differs from r by at least s. 
Moreover, the probability that the super-verifier will pick this Ci and Q is at least 
l/#m = l/poly(n). D 



4 Consistency is QMA-hard 

Theorem 3 Consistency is QMA-hard, via a poly-time randomized Turing reduction 
from Local Hamiltonian. Furthermore, the reduction uses the same value of k for both 
problems, so we get that Consistency with k = 2 is QMA-hard. 

We begin by discussing the basic idea of the proof, and the comphcations that arise. We 
then describe the actual reduction from Local Hamiltonian to Consistency, and finally 
we deal with some issues of numerical precision. 

4.1 The Basic Idea 

Say we are given a local Hamiltonian H = Hi + • • • + Hm, where Hi acts on the subset 
Ci. Consider the following convex program: 

Let p be any 2" x 2" complex matrix. 
Find some p that minimizes tr (Hp), 
such that p ^ and tr{p) = 1. 

It is easy to see that H has an eigenvalue < 7 if and only if the convex program achieves 
ti (Hp) < 7 for some p. (Note that, although the convex program allows mixed states 
p, the optimal p can always be chosen to be a pure state.) Unfortunately, this convex 
program has 4" variables, which makes it unwieldy. 

We now construct another convex program, which is equivalent to the previous one, 
but has only a polynomial number of variables: 

Let pi, . . . , pm be complex matrices, where pi has size 2' *' x 2' *L 
(We interpret each pi as the reduced density matrix for the subset Ci.) 
Find some pi, . . . , pm that minimize ti{Hipi) + • • • + tr{HmPm), 
such that each pi satisfies pi ^ and tr(pj) = 1, 
and pi, ... , Pm are consistent. 

Note that the set of feasible solutions is indeed convex: if pi, . . . , pm are consistent, and 
p[, . . . , p'^ are consistent, then any convex combination p'- = qpi + {^ — q)Pi {i = I, . . . ,m) 
is also consistent. 

Observe that the optimal value of this convex program is equal to the optimal value 
of the previous convex program; this is because, ii pi, . . . , pm are consistent with some 
n-qubit state a, then ti^Ha) = iv{Hipi) + • • • + iv{HmPm)- Also, note that the number 
of variables in this convex program is ^^^^ 4' 'I < 4 ttt, < poly(n). 

This convex program has a "consistency" constraint, which we do not know how to 
evaluate. But if we have an oracle for the Consistency problem, then we can solve this 
convex program, using the algorithm of Bertsimas and Vempala. To make this work, we 
will have to find a suitable representation for the set of feasible solutions, 

K = {{pi, . . . , Pm) which are consistent}. 

Also, we will have to address some questions about the accuracy of the Consistency 
oracle, i.e., how well does it approximate K, and how does this affect the Bertsimas- 
Vempala algorithm. 



We could represent each element (pi, . . . , pm) G i^ by writing down the matrix entries 
for the Pi] then we could view i^ as a subset of C^, where d = Xli^i ^ • -^^^ ^^^^ 
straightforward approach runs into some trouble. Observe that the matrix entries must 
satisfy some algebraic constraints: each pi must be Hermitian, (pj)^ = Pi] and pi and 
Pj must agree on their intersection Q n Cj, that is, trc^-(c,nCj)(Pj) = ^'^Cj-(C,nCj){Pj)- 
Because of these constraints, the set K actually lies in a lower-dimensional subspace of 
C^. We would need to characterize this subspace, before we can apply the Bertsimas- 
Vempala algorithm. We can avoid this problem by switching to a different representation 
for the set K. 

4.2 The Actual Reduction 

We will represent each element of K using the expectation values of the "local" Pauli 
matrices on the subsets Ci, . . . , Cm- These local Pauli matrices form a basis for the space 
of all local Hamiltonians (acting on the subsets Ci). For an n-qubit state o", knowing 
the expectation values of these Pauli matrices is equivalent to knowing the projection of 
a onto this subspace; and this is equivalent to knowing the local density matrices of a. 

First, some notation. Let P be an n-qubit Pauli matrix, and define the "support" 
of P be the set of qubits on which P acts nontrivially; that is, supp(P) = {i | Pj / /}. 
Also, for any subset of qubits C, define the "restriction" of P to C, P\C = (S>iec ^i- 

Define S to be the set of "local" Pauli matrices: 

S = {P e P^" I supp(P) C Ci for some i] - {/}, 

where we excluded the identity matrix / because its expectation value is always 1. Also 
let d=\S\, and note that d < A^m — 1 < poly(n). 

For each P G S, let ap be the corresponding expectation value; and let (ap)pG5 
denote the collection of these ap. Also, let a/ = 1. We define the set K' to be 

K' = {{ap)p(zs which are consistent}, 

where we say the ap are "consistent" if there exists an n-qubit state a such that for all 
P G 5, ap = tr(P(T). Note that K' is a subset of M'^. Also, clearly K' is convex. 

Lemma 4 There is a linear bijection between K and K' . 

Proof: Given some (pi, . . . ,Pm) G K, we can construct {ap)p<zs G K' as follows: 

For each P ^ S: We know that supp(P) C d for some i. So we can write P 
in the form P = (P|Cj) (8) /. Then we set ap = tr((P|Ci)pi). 

If the Pi are consistent with some n-qubit state o", then the ap are also consistent 
with (t; to see this, write ap = tr((P|Cj)/9i) = tr(P(T). (Note that in the case where 
supp(P) C Cj n Cj, it makes no difference whether we pick i or j in the above procedure, 
because pi and pj yield the same reduced density matrix on Cj n Cj.) 

Going in the opposite direction, given some {ap)p^s G K' , we can construct (pi, . . . , pm) G 
K as follows: 



For i = 1, . . . ,m: We construct pi by using the ap for all P with supp(P) C 
Cj. Note that we can write P in the form P = {P\Ci) /. We set 

Pi = ^icT] Yl ap{P\Ci). 

P:supp{P)CCi 

If the ap are consistent with some n-qubit state a, then the pi are also consistent with 
a; to see this, write a in terms of the ap, where we now include the expectation values 
ap = tr(Pa) for ah P G P®", 



i E «-^' 
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note that when we trace out the qubits not in Cj, we get that tr^i^ „}._(7. (P) equals 
2"--|Cil(P|C'j) if supp(P) C d, and otherwise; thus we have 

tr{i,...,n}-c,(o") = ^|c^ E «p(-P|Ci) = Pi- 

P:supp(P)CC, 

Finally, observe that these maps (between K and K') are linear, and they are inverses 
of each other. D 

So we can restate our convex program, using the set K' instead: 

Let ap (for P G 5) be real numbers. 
Find some ap that minimize 



ra 

E^ E aptr(/7,(P|Q)), 
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i=l P : supp(P)CCi 

such that (ap)pg5 G K' (i.e., the ap are consistent). 

Lemma 5 The optimal value of this convex program is equal to the sm^allest eigenvalue 
of the local Ham,iltonian H = Hi + • • • + Hm ■ 

Proof: This follows from the remarks in the previous section, and Lemma HI D 

Next, we prove some bounds on the geometry of the set K' C M*^. 

Lemma 6 K' is contained in a ball of radius R = \fd centered at the origin. 

Proof: Suppose {ap)p(zs G K' , and say it is consistent with some state a. Since ap = 
tr(P(T), it follows that —1 < ap < 1, which implies the result. D 

Lemma 7 The ball of radius r = 1/Va around the origin is contained in K' . 



Proof: Let {ap)p£s be any vector in R of length at most l/vd- By the Cauchy-Schwartz 
inequality, J2p£S \'^p\ — ^' ^^* P ~ "^PeS \'^p\- Now define a = (l/2")(/ + ^pg^ apP). 
This is a legal density matrix, because it can be written as 



a = i^((l-p)/+^(|ap|/ + apP) 



Pes 

. I sr^ , I + sign(ap)P 
= (! — ») h > \ap\ , 

which is (with probability 1 — p) the fully mixed state, and (with probability |ap|, for 
P £ S) the mixture of all eigenstates of P with eigenvalue sign(ap). Furthermore, the 
ap are consistent with a; thus we conclude that {ap)p^s £ K' . ID 

4.3 Numerical Precision 

First, we will show that the Consistency oracle gives a good approximation to the set 
K'. We start by defining a new problem. Consistency', using the expectation values ap 
of the local Pauli matrices P £ S (similar to the definition of K'): 

As in the original Consistency problem, we have an n-qubit system, and 
subsets Ci, . . . , Cm, with \Ci\ < k. But instead of the local density matrices 
pi, . . . , pm, we are given real numbers ap for all P £ S. Each ap is specified 
with poly(n) bits of precision. 

In addition, we are given a real number /?' (specified with poly(n) bits of 
precision) such that /?' > l/poly(n). 

The problem is to distinguish between the following two cases: 

• There exists an n-qubit state a such that, for all P £ S, tv{Pa) = ap. 
In this case, output "YES." 

• For all n-qubit states a, (X]pg5(tr(Pcr) — ap)^) > P'- In this case, 
output "NO." 

Lemma 8 There is a poly-time mapping reduction from Consistenci/ to Consistency. 

Proof: The reduction is as follows: Use the ap to construct pi,. . . , Pm as described in 
Lemma m Set /? = /?'/\/d, where d = \S\. 

Clearly, a "YES" instance of Consistency' maps to a "YES" instance of Consistency. 
Now suppose we have a "NO" instance of Consistency'. Then for all a there is some 
P £ S such that | tT{Pa) — ap\ > P' /\/d. We know that supp(P) C d for some i, so we 
can write P = P ®I where P acts on the subset Cj. Also, let a = trji ,,, ,i}_c^((t). Then 
we have 

|tr(PCT)-tr(Ppi)| >/?'/\/d. 

We will use P to construct a POVM measurement that distinguishes between a and 
Pi. Since the eigenvalues of P are all ±1, we can write P = Hi — 112, where IIi and 112 
are projectors on orthogonal subspaces, and Hi -|- 112 = I- Thus we can use {Hi, 112} as 



a POVM. For the state a, let Sj be the probabihty of measuring j (for j = 1, 2); and for 
the state pi, let rj be the probability of measuring j (for j = 1,2). 
Then we have 

|tr(Pcr) - tiiPpi)\ = \{si- S2) - in - ra)! = 2|si - ri|. 

Observe that the ii distance between s and r is ||s — r||i = |si — ri| + |s2 — ^2] = 2|si — ri|. 
Also, ||s — r||i is a lower bound for the Li (matrix) distance between a and pi. So we 
have 

\\d- - PiWi > \\s - r\\i > f3'/Vd = (3. 

Thus we have a "NO" instance of Consistency. D 

Next, we will show that the Bertsimas-Vempala algorithm succeeds in solving our 
convex program, even when the oracle for the set K' is slightly inaccurate. (The shallow- 
cut ellipsoid method would also work, see [12] for details.) We will make some general 
remarks about the algorithm, and then show that it works for our specific problem. 

For our purposes, we only need to solve a simpler problem, deciding the feasibility 
of a convex program: 

As before, let K he, a, convex set, and let / be a convex function. 
Given some t G M, does there exist a point x £ K such that f{x) < t? 

The Bertsimas-Vempala algorithm is built around a subroutine that solves the fea- 
sibility problem [9]. The basic idea is as follows: 

Let P be the set K. 

Randomly sample some points from P, and compute an approximate centroid 

of P; call this point z. 
If f{z) < t, stop and return true. 

Compute Vf{z), and use this to cut out a portion of the set Pu 
Repeat the procedure starting from line 2. If P gets too small, stop and 

return false. 

The critical step is to sample random points from the set P. (Note that P is convex, 
and we have a membership oracle for P.) One way is to do a random walk known as 
the "ball walk": 

Pick a point y uniformly at random in the ball of radius 5 centered at the 
current position x. If y G P, then move to y, otherwise stay at x. Repeat. 

The points where the membership oracle makes mistakes all lie close to the boundary 
of P; call this the "boundary layer" Pi,. Intuitively, if the boundary layer is thin, it should 
not have much effect on the random walk. Using an argument by Lovasz and Simonovits 
|16j . one can prove (omitting some details): 



^Specifically, we can deduce a hyperplane that separates z from the set {a; | f{x) < t}. Then we talte 
the intersection of P with the half-space that does not contain z. 
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Lemma 9 For any polynomial p, there exists a polynomial q such that, if we run the 
n-dimensional ball walk for at most p{n) steps, and vol{Ph) / vol{P) < l/q{n), then with 
probability 2/3 we will never enter the region Pf,. 

So, if we can show that the boundary layer is small compared to the total volume of 
P, then our algorithm will work fine. (As long as the random walk does not enter 
the boundary layer, the algorithm will perform exactly as if it had access to a perfect 
membership oracle.) 

Finally, there may still be errors due to finite numerical precision — using n bits of 
precision, we have errors of size 2^". This will not be a problem for us, since we only 
need accuracy of l/poly(n). 

We are now ready to prove that the Bertsimas-Vempala algorithm works for our 
specific problem: 

Proof of Theorem O Given an instance of Local Hamiltonian, use Lemma [5] to express 
it as a convex program over the set K' C M'^, where d < poly(n). Let / denote the 
objective function, and set t = {a + h)/2. By LemmaEl we can assume we have an oracle 
for Consistency', which approximates the set K' with error /?', for any (3' > l/poly(n). 
Then use the Bertsimas-Vempala algorithm to solve the following problem: does there 
exist a solution a ^ K' such that f{a) < t? 
Recall that the objective function 



m ^ 



21' 

i=l P : supp{P)CCi 

is linear. We claim that its derivatives in all directions are at most 4^m < poly(n). To 
see this, note that there are at most A^m terms in the sum, and for each term, we have 

I tr(F,P)| < tv{\H,P\) < ||F,||2||P||2 < 2l^»l, 

using the Cauchy-Schwartz inequality for the L2 matrix norm [T7] , and the fact that the 
eigenvalues of Hi lie in the interval [0, 1], while the eigenvalues of P are ±1. 

Now suppose we have a "YES" instance of Local Hamiltonian. Then there exists 
some a* S K' such that f{a*) < a. We claim that the set {a €z K' \ f{a) < t} contains a 
ball of radius 6 > 1/ poly(n). To see this, let a* be the density matrix which corresponds 
to a*. Perturb a* by mixing it with the state //2", then add a small contribution of 
each of the Pauli matrices P G 5. This generates a ball contained in K'. Moreover, 
this ball can have radius 6 > l/poly(n) and still satisfy the condition f{a) < t; this is 
because / does not vary too quickly, and there is a gap between a and t. 

So, in the Bertsimas-Vempala algorithm, the set P always contains a ball of radius 
6. Now set the error threshold for the membership oracle to be /?' < S/d"^, for some 
constant c (to be specified later). We will show that the boundary layer P^ is small 
compared to the total volume of P. Define P^ to be the set P expanded by an amount 
/?', that is, P+ = P + 13' B, where B is the unit ball. We have that 

P+ c P + {P'/5)P = (1 + I3'/5)P, 
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where the equahty holds because P is convex. This imphes that 

vol(P+) < {l + P'/6)'^vol{P) < {l + 2/d^-^)vol{P). 

So we can conclude that vol(Pb) < vol(P+) - vol(P) < {2/d''~^) vol(P). 

Finally, we choose the constant c. We will show, later, that the Bertsimas-Venipala 
algorithm runs in time poly(n). (This upper bound will not depend on the oracle error 
P' , hence it will not depend on c.) We then choose c sufficiently large so that Lemma 
applies. Therefore, the algorithm will work correctly in this case. 

Now suppose we have a "NO" instance of Local Hamiltonian. Then for all a £ K' , 
f{a) > b. In addition, there is some 6 > l/poly(n) such that, for all a within distance 5 
of K' , f{a) > t] this is because / does not vary too quickly, and there is a gap between 
b and t. 

Set the error threshold for the membership oracle to be /?' < 6. Then the set 
{a £ K' \ f{a) < t} is empty, even when the membership oracle makes mistakes. So the 
Bertsimas-Vempala algorithm will work correctly in this case. 

Finally, we claim that the Bertsimas-Vempala algorithm runs in time polynomial 
in n. This follows from Theorem [1] and Lemmas [6] and [71 note that L = \og{R/r) = 
log(poly(n)) = O(logn). D 

5 Discussion 

Consistency of local density matrices is an interesting problem that gives some new 
insight into the class QMA. The reduction from Local Hamiltonian is nontrivial, and in 
that sense. Consistency seems to be an easier problem to deal with. One direction for 
future work is to try to find additional QMA-complete problems by giving reductions 
from Consistency (rather than from Local Hamiltonian). 

Another question is whether Consistency remains QMA-hard under mapping reduc- 
tions. We mention that we can build zero- knowledge proof systems for Consistency 
[18] . using techniques developed by Watrous [19]. If we could show that Consistency is 
QMA-hard under mapping reductions, then we could get zero- knowledge proof systems 
for any language in QMA. 
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